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Abstract 

This paper discusses the possibihty to observe a signal from MSSM Higgs boson 
decays into final states containing four b-jets. Two specific channels are considered: 
bbH and bbA production with H, A — bb, for large values of niH, niA and tan/3, 
and H — > hh ^ bbbb decays for 150 GeV < nin < 2mt and for low values of tan (5. 
Both channels are difficult to extract because of the very large reducible and irre- 
ducible QCD backgrounds. Even with an ultimate integrated luminosity, expected 
per LHC experiment, of 3 • 10^ pb~^, the region of the MSSM parameter space cov- 
ered by these channels does not extend the reach beyond that accessible to other 
channels that were studied in the past. Nevertheless, their observation would help 
in constraining the couplings and branching ratios of the MSSM Higgs bosons. 
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1 Introduction 



One of the most attractive theoretical extensions of the Standard Model (SM) is the 
Minimal Supersymmetric Standard Model (MSSM)[|I], 0. Recently, the interest in 
this model has grown for several reasons. It is the simplest version of low-energy 
super symmetry, with the minimal gauge group and particle content compatible with 
phenomenology. It preserves all SM predictions at low energies, and is therefore not 
contradicted by any of the precision measurements at LEP [^. It also has a high 
level of predictivity in the Higgs sector. At tree level, all Higgs-boson masses and 
couplings can be expressed in terms of only two parameters, usually chosen to be the 
ratio of the vacuum expectation values of the two Higgs doublets, tan jS, and the mass 
of the A boson, mA. Radiative corrections |Q introduce a dependence of the masses 
and couplings on the top-quark mass, the squark masses and the mixing parameters 
in the stop-sbottom sector. All couplings of the MSSM Higgs bosons to fermions 
and bosons can be obtained from the SM couplings by multiplying the latter with 
the appropriate MSSM correction factors, which depend on the parameters listed 
above. 

The MSSM predictions concerning the Higgs-boson masses have important phe- 
nomenological consequences. One of the Higgs bosons, the scalar h, is relatively 
light, with a maximum allowed mass value varying from 105 to 153 GeV, depending 
on the MSSM parameters and the top-quark mass. The masses of the other Higgs 
bosons are expected to be larger and nearly degenerate for mA > 200 GeV. 

There are a variety of production and decay channels through which the MSSM 
Higgs bosons can be observed at LHC Some of these are similar to the SM, such 
as the H, h —>■ 7 7 or H ^ 4£ channels. Others are allowed in the SM but have 
large enough rates only in the MSSM case, e.g. H, A — rr, /i/i, and tt. There are 
also channels characteristic of the MSSM, such as H^ tu, A Zh, H — ^ hh, as 
well as the production of h bosons in squark and gluino decays. 

The work described in p presents an extended and complete study of most 
of these channels, leading to realistic predictions for the LHC discovery potential 
in the MSSM Higgs sector. Figures |l| and ^ show the final results of ||^, namely 
the 5cr-discovery contour curves, for mt = 175 GeV and for integrated luminosities 
of 3 ■ 10'* pb~* and 3 ■ 10^ pb~^ per LHC experiment, respectively, which represent 
the integrated luminosities expected after three and ten years of operation. With the 
initial, modest integrated luminosity of 3 ■ 10* pb~*, a large fraction of the parameter 
space is already covered. With the ultimate integrated luminosity of 3 ■ 10^ pb~*, the 
LHC discovery potential maps the complete parameter space. For the vast majority 
of cases, the experiment would be able to distinguish between the SM and MSSM 
cases. 

The MSSM Higgs sector is quite challenging experimentally for the LHC 0, 
since, most often, the signal-to-background ratios are much smaller than unity and 
the detector resolution in several of the accessible channels is far from optimal. This 
sets stringent requirements on the performance in terms of energy and momentum 
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resolution and of particle identification; however, on the other hand, the variety of 
channels makes this sector an excellent benchmark to evaluate the discovery poten- 
tial of the detector. A very good performance of the electromagnetic calorimeter, 
excellent b-tagging capabilities, good E™^^-resolution and r-identification are all 
crucial ingredients to fully explore the MSSM Higgs sector. 

Channels with b-jets identified in the final state were considered so far only if 
associated with hard leptons or photons: Wh with h — > bb, H/A ^ tt — > £t/jjbb, 
A — ^ Zh ^ iihh and H — >■ hh — ^ bb77. All of these channels were found to be ob- 
servable only for low values of tan/3. The possibility to improve the sensitivity 
by using identified b-jet spectators was also recently considered in the case of 
bbA,bbH production with H/A tt decay. In all these cases, however, isolated 
leptons or photons provide a straightforward experimental trigger, and the purely 
hadronic background from QCD multi-jet production is not a source for concern. 

Final states containing more than two b-jets have not been explored so far in 
these searches. They nevertheless deserve some attention, since they correspond to 
the dominant production and decay modes of the MSSM Higgs bosons over a large 
fraction of the parameter space. Also, it has been reported in recent theoretical 
papers that they should be very promising for discovery in the MSSM Higgs 
sector. 

Two possible channels will be discussed in this paperQ, both leading to final 
states containing four b-jets but no lepton nor photon for the experimental trigger. 

• The first channel arises from bbH and bbA production, which is strongly en- 
hanced for large values of tan/3, since the production rates grow like tan^/3. 
These production processes were studied so far through the tt and /i/i de- 
cay modes over the mass range from 100 to 500 GeV. The tt channel re- 
quires excellent r-identification to suppress the huge background of hadronic 
jets from various sources, but also excellent E^'^'^-resolution for the recon- 
struction of the TT invariant mass. The region of parameter space accessi- 
ble to the TT channel, which has a branching ratio of ~ 10%, is superior 
to that of the /i/i channel; for mA = 300 GeV, the Sa-discovery contour ex- 
tends down to tan/3 ~ 10 for an integrated luminosity of 3 ■ lO'^ pb~^, and 
down to tan/3 ~ 7 for ATLAS+CMS combined with an integrated luminos- 
ity of 3 ■ 10^ pb~^ per experiment. For larger masses, the sensitivity degrades 
rapidly, due to the decreasing production rates and the degradation of the 
TT mass resolution. For mA = 500 GeV, the accessible region reaches only 
down to tan /3 ~ 25 (resp. ~ 18) for an integrated luminosity of 3 ■ 10"^ pb~^ 
(resp. 3 • 10^ pb^^). The rates expected for the bb decay mode would be much 
higher, since the branching ratio is about 90%. However, the observation of 
final states containing multiple b-jets is very challenging: it requires excellent 
b-tagging performance and very efficient jet reconstruction. 

^ The present paper is a slightly modified version of 
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• The second channel is the H — hh decay, which is dominant for low values 
of tan/3 and for ran < 2mt. The observation of this channel would be very 
interesting, since it would correspond to the simultaneous discovery of two 
MSSM Higgs bosons. The final state H — >• hh — bb77 has been studied so far. 
It provides a straightforward experimental trigger and offers good kinematical 
constraints and mass resolution for the reconstruction of mH. The expected 
rates are, however, very low, even when requiring only one identified b-jet in 
the final state. The sensitivity to this channel decreases as tan /? increases, and 
the accessible region of parameter space extends up to tan/3 ~ 2.5 for an in- 
tegrated luminosity of 3 • 10"^ pb~^ and up to tan /3 ~ 4.5 for ATLAS+CMS 
combined with an integrated luminosity of 3 ■ 10^ pb~^ per experiment. The 
H ^ hh ^ bbbb channel would yield much higher rates, by a factor of al- 
most 1000; however, the 4-jet final state will be difficult to trigger on and to 
extract from the huge QCD background. 

The main difficulty and uncertainties in the analysis presented here arise from the 
evaluation of the QCD mult i-jet/ mult i-b- jet backgrounds. There are many sources 
of uncertainties in the predictions for QCD multi-jet production: the parton density 
functions, the choice of scale for the strong coupling constant, the calculations of the 
exact and complete matrix elements and the higher-order corrections. At tree level, 
matrix elements for all parton processes with at most five partons in the final state. 



have been calculated and embedded in a Monte Carlo generator called NJETS |T0 



However, the parton-fiavour information is not accessible in this generator. For 
the studies presented here, the PYTHIA Monte Carlo generator |Tl|, with its QCD 
parton-shower modelling and gluon splitting into heavy quarks, has been used, since 
the flavour content of the jets in the final state is of crucial importance for the evalu- 
ation of many of the background processes. Only leading-order terms are controlled 
in this approach. Predictions from both generators for events with at least 2, 3, 4 
reconstructed high px jets (p^^* > 30, 40, 50 GeV) were compared and results from 
PYTHIA were found to be systematically lower by a factor 2-3. The uncertainties 
on the predictions for high-px multi-jet final states are therefore estimated to be 
as high as at least a factor of 3 and the results presented throughout this paper as 
background rates should be considered as an optimistic. 

The expected detector performance is simulated with the ATLFAST package \T2 . 
This package, used extensively for fast simulation of the ATLAS detector response 
provides reliable estimates of the detector response to hadronic jets. Though the 
parameterisations of the detector performance are ATLAS-specific, the results pre- 
sented here can be considered as representative of what is to be expected at the 
LHC. For those results evaluated in the case of high-luminosity operation, pile-up 
effects are included. The reconstructed jet energies, after kinematical selection, are 
recalibrated, as explained in |12[, to obtain the correct peak position for e.g. the 
m^^ distribution. For the analysis presented here, a rather optimistic b-tagging 
performance is assumed (see for the most recent ATLAS results); for the low- 
luminosity case, an overall b-tagging efficiency of eb = 60% is assumed, with a rejec- 
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tion against c-jets of Rc = 10 and against light-quark and gluon jets of Rj 
For the high-luminosity case, the b-tagging efficiency is degraded to eb 
assuming the same rejection against non-b-jets. 
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Figure 1: For mt = 175 GeV and an integrated luminosity of 3 ■ 10*^ pb~^, expected 
LHC (ATLAS + CMS) 5a-discovery contour curves in the (niA, tanf]) plane for all 
MSSM Higgs-hoson channels studied in 
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Figure 2: For mt = 175 GeV and an integrated luminosity of 3 ■ 10^ pb~^ , expected 
LHC (ATLAS + CMS) 5a -discovery contour curves in the (m\, tan/Jj plane for all 
MSSM Higgs-boson channels studied in j^. 
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Two questions arc of particular interest in the context of these studies: 

• Is it possible, using the bbH/A channel with H/A bb, to extend the region 
in parameter space already covered by the rr channel, particularly to values of niA 
larger than 500 GeV? 

• Is it possible, using the H — > hh — > bbbb channel, to extend the region in 
parameter space already covered by the H — > hh — > bb77 channel? 



2 Observability of H/A — bb for large tan/3 
2.1 Signal events 

For large values of mn and mA, the final state in this channel has a very characteristic 
topology: the two hardest jets in the event come from the H/A — > bb decay, while 
the softer ones come from the associated bb production and from initial/final-state 
radiation. These features have been used for the event selection. At least four 
reconstructed jets are required in the final state, with the three or four jets of highest 
transverse energy tagged as b-jets. The two jets with the largest transverse energy 
are used for the reconstruction of the Higgs-boson mass, m^^g. The selection criteria 
have been optimised for each value of mA and mn, chosen from 300 to 900 GeV, as 
follows: 

• at least four reconstructed jets (ordered in decreasing px) with pt-"^ > 50 GeV 
and pt'''' > 30 GeV; the thresholds on the two hardest jets are raised to: 

• pt^^ > 100 GeV and px^' > 70 GeV (selection S300); 

• pt^i > 200 GeV and pt^" > 100 GeV (selection S500); 

• pt^^ > 250 GeV and pt^" > 150 GeV (selection S700); 

• pt^^ > 300 GeV and pt-"" > 200 GeV (selection S900); 

• the three highest-px jets are required to be tagged as b-jets; an additional 
tagged b-jet is required in two different approaches: 

• in algorithm A, the fourth highest-px jet is chosen; 

• in algorithm B, any tagged b-jet among the additional jets is accepted; 

• the two highest-px jets are used to reconstruct the m^^ distribution; 

• the events are accepted if m^b falls within ±80 to 100 GeV of the considered 
Higgs-boson mass; the mass window was chosen to accept approximately 70% 
of the signal events. 
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The acceptance of the selection 5*500 increases from 10% to 40% as the Higgs- 
boson mass increases from 500 to 900 GeV, as shown in Table ^ The three highest- 
Pt jets are found to be true b-jets only in 6% to 14% of the events, due to initial/final- 
state radiation. This proportion decreases to 1.6% to 3.2% of the events if the fourth 
highest-pT jet is also required to be a true b-jet (algorithm A). For algorithm B, the 
acceptance is improved by about 60%. These rather low acceptances for the signal 
events are unavoidable to reduce the huge QCD backgrounds to an acceptable level, 
as discussed in Section 2.2. It was checked for example that the fact of requiring 
any combination of the reconstructed jets to be identified as b-jets leads to much 
worse signal-to-background ratios and significances. 

Table ^ shows, for an integrated luminosity of 3 ■ 10^ pb~^, the expected rates 
of signal events after applying the selection 5*500 and the b-tagging procedure, for 
different assumptions on the b-tagging performance in addition to the default one 
described above. The signal sample is dominated by events containing true b-jets: 
the contribution from events with at least one mis-identified jet is below 10%, as 
shown by comparing the results for the default jet rejections with those for an almost 
infinite rejection, R = 10^. 




Figure 3: Reconstructed m^s mass after applying the selection 5*500, for the bbH 
signal with H ^ bb decay and for mn = 500 GeV. The degradation in resolution 
expected from pile-up at high luminosity has been included. 
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Table 1: Acceptance of selection 6*500 and fractions of events containing one or more 
true b-jets, for three values o/niA and hih. 



bbH, bbA 


a 


Four 


Jetl 


+ jet2 


+ jet3 


Algorithm 


Algorithm 


with H/A ^ bb 


(pb) 


jets 


= true 


= true 


= true 


A 


B 


Selection S500 






b-jet 


b-jet 


b-jet 






mn = 500 GeV 


3.5 


10.8% 


10.0% 


8.4% 


5.6% 


1.6% 


2.4 % 


niH = 700 GeV 


0.9 


20.4% 


18.0% 


14.8% 


8.0% 


2.0% 


3.2 % 


niH = 900 GeV 


0.2 


39.2% 


36.0% 


26.4% 


13.6% 


3.2% 


5.5 % 



Table 2: Expected rates of reconstructed signal events after applying selection S500 
and b-tagging, for tan f3 = 30 and for an integrated luminosity of 3 ■ 10^ pb~^ (AT- 
LAS). The results are shown for three different assumptions on the b-tagging per- 
formance. 



Selection S500 


Four 


Jetl 


+ jet2 


+ jet3 


Algorithm 


Algorithm 




jets 


tagged 


tagged 


tagged 


A 


B 


niA, niH 




as b-jet 


as b-jet 


as b-jet 










eb = 0.5 


, ec = 10-^ R = 10« 


500 GeV 


111000 


52400 


22000 


7250 


1090 


1670 


700 GeV 


54600 


24700 


9850 


2660 


363 


572 


900 GeV 


23800 


10700 


4030 


1030 


125 


217 






eb = 0.5, ec = 


0.1, R=10^ 


500 GeV 


111000 


52500 


22100 


7400 


1180 


1770 


700 GeV 


54600 


24700 


9940 


2750 


404 


622 


900 GeV 


23800 


10700 


4080 


1070 


142 


238 






eb = 0.6, ec = 


0.1, R=10^ 


500 GeV 


111000 


63000 


31900 


12700 


2420 


3630 


700 GeV 


54600 


29700 


14300 


4730 


824 


1270 


900 GeV 


23800 


12900 


5860 


1830 


288 


486 



The reconstructed mbb peak is rather broad, because final-state radiation and 
hadronisation degrade the mass resolution and also because a significant fraction 
(about 20%) of the bb combinations entering the distribution are incorrect (one or 
both b-jets do not originate from the H — >^ bb decay). Figure ^ shows as an example 
the distribution of the reconstructed mbb ^^i mn = 500 GeV; the acceptance in the 
chosen mass bin (mn ± 80 GeV) is about 70%. 

The quality of the mass reconstruction has been studied separately, using the 
gg — s> H — s> bb process, which does not suffer from the combinatorial backgrounds 
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of the bbH process. The mass resolution obtained and the acceptance in the mass 
window are shown in Table ^ and in Fig. ^ for all events without any selection. The 
large tails in the mass distributions that can be seen in Fig. ^ are almost completely 
removed if one applies the selection criteria on the jet transverse energies, but the 
events in the tails are nevertheless lost. Most of these losses are accounted for in 
the acceptances quoted in Table |I|. Finally, Table | shows the total expected signal 
rates for three values of mn and the appropriate selection cuts, before and after 
applying the chosen m^s mass window (including pile-up at high luminosity), for an 
integrated luminosity of 3 ■ 10^ pb~^. 



Table 3: Mass resolutions expected for the reconstructed m^^ mass, for the 
gg — H — * bb process in the absence and presence of pile-up. No selection cuts 
are applied. 



Higgs mass (GeV) 


No pile-up 


With pile-up 




am (GeV) 


Within ± 2 dm 


(Tn, (GeV) 


Within ± 2 (Tm 


300 


22.7 


62% 


26.4 


60% 


500 


33.8 


57% 


39.7 


55% 


700 


38.3 


52% 


44.0 


50% 


900 


52.6 


53% 


63.5 


52% 



Table 4: Expected numbers of signal events after applying the appropriate selec- 
tion criteria, for an integrated luminosity of 3 ■ 10^ pb~^ (ATLAS). The numbers of 
events accepted within the m^b mass window are shown in brackets. 



bbH, bbA production 


Kinematical 


Three 


Algorithm 


Algorithm 


with H/A ^ bb 


cuts 


b-jets 


A 


B 


mn = 500 GeV 


1.1 ■ 10^ 


7400 


1100 


1800 


Selection S500 




(5000) 


(960) 


(1300) 


mn = 700 GeV 


3.4 • 10^ 


2289 


360 


530 


Selection S700 




(1530) 


(240) 


(350) 


mn = 900 GeV 


1.3 ■ 10^ 


900 


135 


210 


Selection S900 




(524) 


(80) 


(122) 
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Figure 4: Reconstructed 111^5 before any selection cuts, for the gg ^ H ^ bb process 
and for hih = 300, 500, 700 and 900 GeV. The degradation in resolution expected 
from pile-up at high luminosity has been included. 

2.2 Background events 

The background to the bbH/A signal with H/A — > bb decay is completely dominated 
by QCD multi-jet processes with final states containing a variable number of real 
b-jets. For the analysis presented here, large samples of background events were 
generated using PYTHIA, based on the hard-scattering di-jet sub-processes. The 
hard-scattering process was accompanied by initial-state radiation, which provides 
one of the sources of additional jets; gluon splitting and/or final-state radiation also 
provide a source of additional jets, some of which originate from b-quarks. 

The event generation was organised in several px bins with different statistics, 
typically 5 ■ 10^ events to 4 ■ 10^ events per bin. A sufficiently large number of events 
were generated in each bin to obtain statistically significant samples after applying 
the selection criteria for the final states of interest. 

As for the signal, the background events were simulated through the ATL- 
FAST package, with pile-up effects included for the high-luminosity case. The b- 
tagging procedure and the jet energy recalibration were also applied. Table |^ shows 
the summed production cross-sections and the kinematical acceptances for selection 
5*500 (before applying the b-tagging procedure), for the most prominent background 
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sub-processes: bb, qq, gb, gq and gg. The kinematical acceptances are about 0.13% 
for all the di-jet sub-processes, but the initial cross-sections vary over several orders 
of magnitude. Figure ^ shows the contributions from the various pt bins chosen at 
generation to the distribution of m^^^. The relative contributions from the different 
bins vary rapidly depending on the selection criteria, so the use of this algorithm for 
event generation turned out to be rather efficient in this caseQ. As an example, even 
though the production cross-sections in the first two bins are different by several 
orders of magnitude (2 ■ 10^ pb and 3.5 ■ 10^ pb respectively), the contributions 
to the total background from these bins after selection S500 are comparable. 

After selection, but before applying the b-tagging procedure, the inclusive back- 
ground rates are approximately a factor of 10^ to 10^ higher than the signal rates in 
the mass bin of interest. Requiring at least three identified b-jets reduces this factor 
to 10^ to 10^. Requiring in addition a fourth identified b-jet gains another factor of 
about 2. In the sample containing at least three identified b-jets, approximately 40% 
of the background events contain at least three true b-jets, whereas in the sample 
containing at least four identified b-jets this fraction increases to 65%. These frac- 
tions can be estimated directly from the numbers shown in Table ^ by comparing 
the expected rates for an almost infinite rejection against non-b-jets (R = 10^) with 
the default value of this rejection (R = 10^). 

The dominant remaining background arises from gb and gg production with 
gluon splitting into a bb pair. The contribution from direct gg bb production is 
found to be only at the level of ~ 10% of the total background. The contribution 
from tt production is very small (~ 1%). Figure ^ shows the respective contributions 
from the remaining background processes to the reconstructed m^bhiass distribution, 
namely from direct bb, from the summed bb + g b and from the summed all sub- 
processes, after applying the b-tagging procedure with the performance expected at 
high luminosity. 

Table 5: Production cross-sections and kinematical acceptances of the selection cuts 
described in Section 2.1, for various QCD background processes. 



Sub-process 


a (pb) 


Acceptance for 




pf > 50 GeV 


selection S500 


gg, qq gg, qq 


1.32 ■ 10'^ 


0.12% 


gq gq 


8.23 ■ 10^ 


0.16% 


gb gb 


4.40 • 10^ 


0.15% 


gg, qq qq 


4.32 ■ 10^ 


0.15 % 


gg, qq bb 


8.46 ■ 10^ 


0.14% 


tt 


5.05 ■ 102 


4.20% 


Total 


2.22 • 10^ 


0.13% 



■^Many thanks to F. Paige for this suggestion. 
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Figure 5: Distribution of the reconstructed mass for background events after 
applying the selection cuts S500 (top), S700 (middle) and S900 (bottom), for an in- 
tegrated luminosity of 3 ■ 10^ pb^^ (ATLAS). The contributions from the different 
Pt bins used at generation are shown separately: 50 < px^""^ < 200 GeV (solid), 
200 < ptS""" < 350 GeV (dashed), 350 < pt^""" < 500 GeV (dots), and 
Pt^^" > 500 GeV (dot-dashed). 

The background estimates presented above explain why the theoretical estimates 
of PI are much too optimistic. These estimates were based on an evaluation of the 
QCD background in final states containing three b-jets using only the contribution 
from the direct gg — > bb sub-process. As appears clearly in Table |], this sub-process 
contributes not more than 10 to 15% of the total background for final states con- 
taining at least three identified b-jets. This conclusion was found to be rather 
independent of the selection cuts applied to the reconstructed jets, as demonstrated 
by the results of Section 3.1 in Table |12[ 

Finally, Table ^ shows the total expected background for the different selection 
cuts, before and after applying the chosen m^g window, for an integrated luminosity 
of 3 ■ 10^ pb~^ and the ATLAS experiment. 
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Table 6: Expected rates of reconstructed background events after selection S500 as a 
function of the b-tagging performance, for an integrated luminosity of 3 • 10^ pb~^ 
(ATLAS). 



Sub-process 


Four 


Jctl 


+ jet2 


+ jet3 


Algorithm 


Algorithm 




jets 


tagged 


tagged 


tagged 


A 


B 


Selection 5*500 






as b-jet 


as b-jet 


as b-jet 










eb = 0.5, ec = 10-^ R = 10*^ 


gg, qq bb 


3.5 


10^ 


1.4 


•10^ 


3.7 • 10^ 


7.2 • 10^ 


7.3 • 10=^ 


1.5 


• 10^ 


gg, qq qq 


2.0 


10^ 


1.6 


10^ 


5.2 ■ 10^ 


3.7-10^ 


3.9-103 


1.0 


• 10^ 


gb gb 


1.8 


10« 


4.5 


10^ 


4.2 • 10*^ 


2.3 • 10^^ 


2.2 • 10^ 


4.1 


• 10^ 


gq^ gq 


3.9 


109 


3.5 


• 10^ 


5.5- 10^ 


2.2 • 10^ 


1.0 • 10^ 


8.4 


• 102 


gg, qq gg, qq 


4.8 


10^ 


8.2 


• 10^ 


4.6 • 10^ 


1.8 • 10^ 


1.4 • 10^ 


5.2 


• 10^ 


Total 


9.1 


10« 


1.9 


•10« 


1.8-10^ 


5.4 • 10^ 


4.7 • 10^ 


1.2 


• 10^ 








ri, = 0. 


'). r, = 0.1. R = 10 






gg, qq bb 


3.5 


10^ 


1.4 


■10^ 


3.9 ■ 10"^^ 


1.5 ■ lO'^ 


1.0 ■ 10^ 


1.9 


• 10'^ 


gg, qq qq 


2.0 


10^ 


2.0 


■10^ 


5.4- 10^ 


1.2 ■ 10^ 


6.1 ■ lO-'^ 


1.5 


■ 10^ 


gb gb 


1.8 


108 


4.6 


•10^ 


5.1 • 10^ 


3.9 • 10^ 


3.1 • 10^ 


5.5 


• 10^ 


gq^ gq 


3.9 


109 


9.2 


• 10^ 


2.7-10'' 


1.3 • 10^ 


7.0 • 10=* 


1.3 


• 10^ 


gg, qq gg, qq 


4.8 


10^ 


1.5 


•10« 


7.8 • 10^ 


5.0 • 10^ 


2.8 • 10^ 


7.8 


• 10^ 


Total 


9.1 


109 


3.2 


•10« 


2.5 • 10^ 


1.3 • 10^ 


8.2 ■ 10^ 


1.8 


• 10^ 








= 0.6, = 0.1, R = 10 


2 




gg, qq bb 


3.5 


10^ 


1.7 


■10^ 


5.5-106 


2.3 ■ 10^ 


2.0 ■ 10^ 


3.7 


■ 10^ 


gg, qq qq 


2.0 


108 


2.3 


•10^ 


7.7- 10^ 


1.8 • 10^ 


1.1 • 10^ 


2.8 


• 10^ 


gb ^ gb 


1.8 


10« 


5.5 


•10^ 


7.1 • 10^ 


6.3 • 10^ 


6.1 • 10^ 


1.1 


• 10^ 


gq^ gq 


3.9 


109 


9.9 


•10^ 


3.2 • 10^ 


1.9 • 10^ 


1.1 • 10^ 


2.0 


• 10^ 


gg, qq gg, qq 


4.8 


109 


1.7 


10*^ 


1.0 ■ 10^ 


7.5 ■ 10^ 


5.0-10^ 


1.5 


■ 10^ 


Total 


9.1 


109 


3.6 


•10« 


3.4 • lO'' 


2.0 • 10*^ 


1.5 • 10^ 


3.4 


• 10^ 



Table 7: Expected rates of reconstructed background events after applying the ap- 
propriate selection criteria, for an integrated luminosity of 3 ■ 10^ pb~^ (ATLAS). 
The numbers of events accepted within the chosen m-^^^ mass window are shown in 
brackets. 





Kinematical 


Three 


Algorithm 


Algorithm 




cuts 


b-jets 


A 


B 


Selection S500 


9.2 • 109 


1.3 • 10^ 


8.2 • 10^ 


1.8-10^ 






(3.9 ■ 10^) 


(3.3 ■ 10^) 


(6.0 ■ 10^) 


Selection S700 


5.1 • 109 


7.2 ■ 10''' 


4.0 ■ 10^ 


1.2 ■ 10^ 






(1.7 ■ 10'^') 


(1.0- 10^) 


(2.7-104) 


Selection S900 


2.0 ■ 109 


2.9 ■ 10^ 


1.8 ■ 10^ 


3.0-104 






(4.8 • 10^) 


(2.8 • 10^) 


(5.2 - 10^) 
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Figure 6: Distribution of the reconstructed m^^, mass for background events after ap- 
plying the selection cuts S500 (top), S700 (middle) and S900 (bottom) and the b-tagging 
procedure, for an integrated luminosity of 3 • 10^ pb~^ (ATLAS). The separate con- 
tributions from direct bb (hashed), from bb + gb (dashed) and the total spectrum 
(solid) are also shown. 
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2.3 Expected signal significances 

Table || (for an integrated luminosity of 3 ■ 10^ pb~^) and Table ^ (for an integrated 
luminosity of 3 ■ 10^ pb~^) show the final signal and background rates and expected 
significances for ATLAS experiment as a function of the selection procedure adopted, 
for tan (3 = 30 and for four values of the Higgs-boson masses. The selection criteria 
(S300 to S900), as well as the mass window for the reconstructed m^b, were optimised 
separately for each mass (a mass window of ~ ± 2am was used). 

The overall conclusion is that, even for an integrated luminosity of 3 ■ 10^ pb^^, 
the extraction of the signal for bbH/A production with H/A ^ bb decay will be very 
difficult. Although the expected rates and significances are higher for a selection 
requiring only three identified b-jets, it is unlikely that the systematic uncertainties 
on the background shape can be controlled to the required precision in this case, 
where the signal-to-background ratio is lowest (S/B ~ 1%). A selection requiring 
at least four identified b-jets yields a more favourable signal-to-background ratio 
of ~ 3% and also a more favourable ratio (of ~ 67%) of irreducible-to-total back- 
ground0. As mentioned previously, this latter ratio can be easily extracted from the 
numbers given in Table ^ for R = 10^ and R = 10^. 

For a selection requiring four identified b-jets and for an integrated luminosity 
of 3 ■ 10^ pb~^ (ATLAS), a nominal significance larger than 5a could be achieved for 
tan/5 > 29 (mn = 500 GeV) and for tan/3 > 20 (mH = 300 GeV). This coverage 
in the parameter space can be compared with that of the rr channel, which extends 
down to tan/? > 25 (mn = 500 GeV) and to tan/3 > 8 (mn = 300 GeV). For 
an integrated luminosity of 3- 10^ pb~^, the respective lower limits on tan/3 are 
tan/? > 42 and tan /? > 30 for the bb channel, and tan /? > 33 and tan/? > 11 
for the rr channel. 

It should be stressed again that the sensitivities quoted for the bb channel are 
rather on the optimistic side, since the estimates of the QCD background are very 
uncertain (the PYTHIA results could be under-estimated by a factor of 3), and the 
assumptions used for the b-tagging performance are more optimistic than the results 



from recent work for the ATLAS Inner Detector Performance TDR |13|. In addition, 
systematic uncertainties due to the lack of knowledge of the background shape have 
not been taken into account in the significance estimates of Tables H and |. 



^The irreducible background consists of the events containing the required number of identified 
true b-jets. 
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Table 8: Expected signal and background rates inside the uiy^^ mass window, for 
four values of hih and niA, for tan/3 = 30 and for an integrated luminosity 
of 3 ■ 10'' pb~^ (ATLAS). Also shown are the expected signal-to-background ratios 
and signal significances for various selection algorithms (see text). 



bbH/A with H/A ^ bb 


Signal (S) 


Background (B) 


S/B S/VB 




ma = 300 GeV 


Selection S300 










Three b-jets 


3200 


2.1 • 10^ 


1.5% 


7.0 


Algorithm A 


630 


2.4 • 10^ 


2.6% 


4.1 


Al&oiitlnii B 




1075 


1.6 • 10' 


2.3% 


5.0 




mn = 500 GeV 


oeiecLion 0500 










Three b-jets 


855 


5.9 ■ 10^ 


1.4% 


3.5 


Algorithm A 


200 


6.1 • 10^ 


3.3% 


2.6 


Algorithm B 


265 


1.1 • 10^ 


2.4% 


2.5 




mn = 700 GeV 


Selection S700 










Three b-jets 


260 


2.6 • 10^ 


1.0% 


1.6 


Algorithm A 


50 


1.8 • 10^ 


2.8% 


1.2 


Algorithm B 


70 


5.1 • 10-^ 


1.4% 


1.0 




mn = 900 GeV 


Selection S900 










Three b-jets 


90 


7.3 • 10^ 


1.2% 


1.0 


Algorithm A 


16 


5.2 • 10^ 


3.0% 


0.7 


Algorithm B 


25 


9.9 • 10^ 


2.5% 


0.8 
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Table 9: Expected signal and background rates inside the uiy^^ mass window, for 
four values of hih and niA, for tan/3 = 30 and for an integrated luminosity 
of 3 ■ 10^ pb~^ (ATLAS). Also shown are the expected signal-to-background ratios 
and signal significances for various selection algorithms (see text). 



bbH/A with H/A ^ bb 


Signal (S) 


Background (B) 


S/B S/VB 




ma = 300 GeV 


OCicL^ llUii O3QO 










Three b-jets 


18700 


1.4 • 10*^ 


1.3% 


15.8 


Algorithm A 


3080 


1.3 • 10^ 


2.4% 


8.5 


Algorithm B 


5270 


2.4 • 10^ 


2.2% 


10.7 




mn = 500 GeV 


Ot!lt!CtlOIl 0500 










Three b-jets 


5000 


3.9 ■ 10^ 


1.3% 


8.0 


Algorithm A 


960 


3.3 • 10^ 


2.9% 


5.3 


Algorithm B 


1300 


6.0 • 10^ 


2.2% 


5.3 




mn = 700 GeV 


Selection Sooo 










Three b-jets 


1530 


1.7- 10^ 


0.9% 


3.7 


Algorithm A 


240 


1.0 ■ 10^ 


2.4% 


2.4 


Algorithm B 


350 


2.7- lO'^ 


1.3% 


2.1 




mn = 900 GeV 


Selection S900 










Three b-jets 


524 


4.8 • 10^ 


1.1% 


2.4 


Algorithm A 


80 


2.8 • 10-^ 


2.8% 


1.5 


Algorithm B 


122 


5.2 • 10^ 


2.3% 


1.7 
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3 Observability of H — > hh — > bbbb 
for small tan (3 

3.1 Signal events 

Two points in the MSSM parameter space were chosen for a detailed study of this 
channel: mn = 300 GeV with tan j3 = 1.5 and tan (3 = 3.0, corresponding respectively 
to nih = 77.8 GeV and nih = 98 GeV (2-loop calculation). The final state in this 
channel can be fully reconstructed with good mass resolution through the use of the 
constraint, m-^^ = m^. Four identified b-jets are required with two bb combinations 
reconstructed with m^b close to m^. As shown also in |]14|, with the expected ATLAS 



calorimeter performance, a signal acceptance of 70% to 80% can be achieved using 
a mass window, m^b = nih ± 25 GeV, for the h-boson reconstruction and, after 
applying the constraint on m^, a mass window, mbbbb = i^h ± 20 GeV. 

For the signal events, the jet transverse-energy spectrum is rather hard, with 
< pt-*^ > = 100 GeV and < px-'"' > = 36 GeV for jets within the Inner Detector ac- 
ceptance for b-tagging. For the high-luminosity case, at least four jets reconstructed 
with pt > 40 GeV (before energy recalibration) are required, yielding an acceptance 
of ~ 25% for the signal events. Table ^ shows the expected signal rates for two 
different jet-selection algorithms described below: 

• Algorithm A: 

• the four most energetic jets are required to be identified as b-jets; 

• all possible bb combinations are then reconstructed; 

• the best pair of combinations for the reconstruction of mh is chosen by 
minimising = {^hh,i - ^h)^ + (j^hh,2 - ^hY'i 

• both pairs are required to satisfy mbb = nih ± 25 GeV; 

• the event is accepted if mbbbb = ± 20 GeV after applying a constraint 
on mh. 

• Algorithm B: 

• all possible jj combinations inside the Inner Detector acceptance are consid- 
ered; 

• the best pair of combinations for the reconstruction of mh is chosen by 
minimising = (nijj i ~ nih)^ + (mjj ^ — m^Y; 

• all four jets chosen in this way are required to be identified as b-jets; 

• both pairs are required to satisfy mbS = nih ± 25 GeV; 

• the event is accepted if mbbbb = ± 20 GeV after applying a constraint 
on mh. 

A careful comparison of these algorithms with other possible selection methods has 
shown that Algorithm B is close to optimal. 
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Table 10: Numbers of expected H — > hh ^ bbbb signal events as a function of 
the selection algorithm, for an integrated luminosity of 3 ■ 10^ pb~^ (ATLAS). The 
high-luminosity b-tagging performance is assumed (e^ = 50%, ec = 10% and 

R = loo;. 







Four 


Four 


+ m^b within 


+ mn within 


niH = 300 GeV 


a 


reconstructed 


identified 


± 25 GeV 


± 20 GeV 




(pb) 


jets 


b-jets 


of mh 


of mH. 




Algorithm A 


tan/3 = 3.0 


0.76 


5.8 ■ 10^ 


350 


154 


125 


tan/3 = 1.5 


1.73 


1.6 ■ 10^ 


800 


510 


430 




Algorithm B 


tan/3 = 3.0 


0.76 


5.8 • 10^ 


640 


470 


360 


tan/3 = 1.5 


1.73 


1.6 ■ 10^ 


1760 


1560 


1360 



3.2 Background events 

As in Section 2.2, the background events arise dominantly from QCD multi-jet 



production. Table [Tl| shows the expected rates of reconstructed multi-jet back- 
ground events for the various QCD sub-processes and for an integrated luminosity 
of 3 ■ 10^ pb~^. This Table also shows the contributions from events containing one 
or more true b-jets to each of the sub-processes; for the dominant gg sub-process, 
only 0.01% of the events contain four true b-jets in the final state, whereas this 
factor increases to 0.8% for the direct bb sub-process. Nevertheless, since the pro- 
duction cross-sections differ by two orders of magnitude, the direct bb production 
process contributes only ~ 10% of the total QCD background after having required 



four tagged b-jets, as shown in Table 0. The dominant background is from the 



gb sub-process, and the irreducible background from four true b-jets is shown to 
amount to ~ 70% of the total multi-jet background. This can be deduced from a 
comparison of the expected background rates in Table |T2] for R = 10^ (almost infi- 
nite rejection of non- b-jets) and for R = 10^ (default b-tagging performance). This 
is also illustrated in Fig. |^, which displays the mass distributions expected for the 
candidate h bb decays for the three dominant sub-processes and for the events 
accepted by algorithm B. Finally, Table |T3| shows the expected rates of background 
events at each step of the selection procedures described as Algorithms A and B. 
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Table 11: Expected rates of QCD multi-jet background events for different sub- 
processes and for an integrated luminosity of 3 • 10^ pb"^ (ATLAS). Also shown are 
the contributions to the total background from events containing one or more true 
b-jets. 



Sub-process 


Four jets 


Jetl 


+ jet2 


+ jets 


+ jet4 




with pt > 40 GeV 


= true 


= true 


= true 


= true 




and 77 < 2.5 


b-jet 


b-jet 


b-jet 


b-jet 


gg, qq bb 


2.0 • 10^ 


70% 


30% 


2% 


0.8% 


gb gb 


1.2 • 10^ 


45% 


10% 


1% 


0.3% 


gg, qq gg, qq 


1.7-10^° 


3% 


0.4% 


0.04% 


0.01% 



Table 12: Expected rates of reconstructed background events after jet- selection cuts as 
a function of the b-tagging performance, for an integrated luminosity of 3 • 10^ pb"^ 
(ATLAS). 



Sub-process 


Four jets 


Jetl 


+ jet2 


+ jet3 


+jet4 




with Pt > 40 GeV 


tagged 


tagged 


tagged 


tagged 




and \r}\ < 2.5 


as b-jet 


as b-jet 


as b-jet 


as b-jet 




n, = 


0.5. r, = 


IQ-". R 


= 10'' 


gg, qq bb 


2.0 • 10^ 


6.9 • 10^ 


1.5-10' 


5.0 - 10^ 


9.7 - 10^ 


gb gb 


1.2 • 10^ 


2.7- 10*^ 


2.9 - 10' 


1.5 - 10^ 


2.1 - 10^ 


gg, qq gg, qq 


1.7- 10^° 


2.6 ■ 10*^ 


1.6 - 10' 


8.6 - 10^ 


1.2 - 10^ 


Total 


1.8 ■ 10^° 


6.0 ■ 10*^ 


6.0-10' 


2.9 - 10^ 


4.3 - 10^ 




Cb = 


0.5, ec = 


: 0.1, R 


= 10^ 


gg, qq bb 


2.0 • 10« 


7.0 • 10' 


1.6 • 10' 


9.5 - 10^ 


1.2 - 10^ 


gb^gb 


1.2 • lO'' 


2.8 • 10^ 


3.4 - 10' 


2.6 - 10^ 


2.9 - 10^ 


gg, qq gg, qq 


1.7-10^° 


5.2 ■ 10*^ 


2.8 - 10' 


2.1 - 10^ 


1.9 - 10^ 


Total 


1.8 • 10^° 


8.7- 10^ 


7.8- 10' 


5.6 • 10^ 


6.0 - 10^ 




£6 = 


0.6, ee = 


0.1, R 


= 10^ 


gg, qq bb 


2.0 • 10^ 


8.4-10' 


2.2 - 10' 


1.5 - 10^ 


2.4-10^ 


gb gb 


1.2 • 10^ 


3.3- 10^ 


4.8-10' 


4.2 - 10^ 


5.6-10^ 


gg, qq gg, qq 


1.7-10^° 


5.7- 10« 


3.6-10' 


3.2 - 10^ 


3.5 - 10^ 


Total 


1.8 • 10^° 


9.8 • 10^ 


1.1 - 10^ 


8.9 - 10^ 


1.1 - 10^ 
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Figure 7: Distribution of the reconstructed m^b for background events after applying 
selection B, shown separately for the various sub-processes: bb (top), gb (middle) 
and gg + gq + qq (bottom). Also shown separately are the contributions from events 
containing four true b-jets (solid), three true b-jets (dashed), two true b-jets (dots) 
and one or less true b-jets (dot-dashed). 

Table 13: Expected rates of reconstructed background events as a function of the 
selection cuts, for algorithms A and B and for an integrated luminosity of 3 ■ 10^ pb~^ 
(ATLAS). The default high luminosity b-tagging performance is assumed. 



Sub-process 


a (mb) 
Pt > 10 GeV 


Four jets 
with Pt > 40 GeV 
and |?7| < 2.5 


Four jets 
identified 
as b-jets 


nih within 
± 25 GeV 
of nominal 


mn within 
± 20 GeV 
of nominal 




Algoiitlnii A 


bb 

gb 

gg, gq, qq 


2.1 • 10-2 
7.1 ■ 10-2 
5.9 


2.0 ■ 10^ 
1.2 ■ 10^ 
1.7- 10^° 


1.2 ■ 10^ 

2.8 ■ 10^ 

1.9 ■ 10^ 


2000 

13000 

2000 


200 

1900 

500 


Total 




1.8 ■ 10^° 


5.9-105 


17000 


2600 




Algorithm E 


5 


bb 
gb 

gg, gq, qq 


2.1 • 10"^ 
7.1 ■ 10-2 
5.9 


2.0-10** 

1.2 ■ 10^ 
1.7-101° 


1.4 • 10^ 

2.9 ■ 105 
1.9 ■ 105 


3000 

24000 
14000 


500 

2500 
1000 


Total 




1.8- 10^° 


6.2 • 105 


41000 


4000 
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3.3 Significance 

Table shows, for both algorithms A and B, the expected signal and background 
rates for nin = 300 GeV, for tan/3 =1.5 and 3.0, and for an integrated luminosity 
of 3 • 10^ pb~^. Algorithm B yields a significantly better signal-to-background ratio 
and thereby sensitivity, since b-jets from H — > hh — > bbbb decays are not very often 
the hardest jets in the event. 

The m^b and m^bbb mass distributions for both the signal and background events 
selected by algorithm B are shown in Fig. ^ and Figure § before the relevant mass 
cuts are applied. Fig. |TU| shows the distributions of p!^'', p!^'^'^^, p^ and p^ before 



energy recalibration for the accepted signal and background events. After applying 
the selection cuts of algorithm B, these distributions do not differ sufficiently between 
signal and the background to justify further optimisation of the selection cuts. 

A sensitivity similar to that expected for the H hh — > bb77 channel (see Fig. ^) 
appears to be also achievable in this channel, since rate is larger, but this would 
require triggering on four-jet events at high luminosity with p-!^ > 40 GeV and 
|?7| < 2.5. The expected trigger rates are very uncertain, in the range between 1.6 
and 4.5 kHz for PYTHIA and NJETS respectively (before energy recalibration), and 
also too high to be accepted without further cuts at level-2. The possibility of using 
b-tagging to reduce the rate will be investigated and compared with other simpler 
possibilities involving somewhat tighter cuts. 

This channel has also been studied in the low-luminosity case to evaluate the 
possible gain in sensitivity that could be achieved by lowering the jet ET-threshold. 
The results of this study are shown, for algorithm B, for tan (3 = 3.0 and for an inte- 
grated luminosity of 3 ■ 10'^ pb~^, in Table |l^ as a function of the jet Ex-threshold, 
which was varied from its minimum realistic value (before recalibration) of 15 GeV 
to the high-luminosity value of 40 GeV. The background rates decrease much faster 
than the signal rates and the sensitivity improves as the jet ii^T-threshold is raised. 
For a threshold of 20 GeV, Fig. |ll] shows the distributions of p^^, pJjJ''^^, pJ^^ and 



Pt* for the accepted signal and background events: requiring a cut on p^, typically 
> 60 GeV, would obviously improve the signal sensitivity. As an illustration. 
Table ^ also shows the resulting expected signal and background rates as well as the 
signal-to-background ratios and signal significances, for events with p!^^ > 60 GeV 
and with p^' > 80 GeV, as a function of the jet Ex-threshold. 

Even with optimised cuts, the sensitivity at low luminosity is weaker than that 
for the H — >■ hh — >■ bb77 channel 0, and the trigger requirements become very de- 
manding, since an Ex-threshold of about 20 GeV would be desirable. 
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Table 14: Form^ — 300 GeV and two values of tan (3, expected signal and background 
rates, signal-to-background ratios and signal significances, for the two selection algo- 
rithms A and B described in the text and for an integrated luminosity of 3 ■ 10^ pb~^ 
(ATLAS). 



H ^ hh ^ bbbb 


a(pb) 


S 


B 


S/B 


s/Vb 




Algorithm A 


tan/3 = 3.0 


0.76 


125 


2600 


4.8 % 


2.4 


tan /3 = 1.5 


1.73 


430 


2600 


16.5 % 


8.4 




Algorithm B 


tan/3 = 3.0 


0.76 


360 


4000 


9.0 % 


5.7 


tan /3 = 1.5 


1.73 


1360 


4000 


34.0 % 


21.5 



Table 15: For mn — 300 GeV and tan /3 — 3.0, expected signal and background rates, 
signal-to-background ratios and signal significances as a function of the chosen jet 
'Ei^ -threshold, for an integrated luminosity of 3 ■ 10^ pb~^ (ATLAS). 



H ^ hh ^ bbbb 


S 


B 


S/B 


s/Vb 


Jet Et threshold 


Algorithm B 


15 GeV 


275 


64 000 


0.4 % 




1.1 


20 GeV 


231 


48 000 


0.5 % 




1.1 


30 GeV 


132 


10 000 


1.3 % 




1.3 


40 GeV 


50 


800 


6.2 % 




1.8 




+p!jf' > 60 GeV and 


> 80 GeV 


15 GeV 


160 


3 900 


4.1 % 




2.6 


20 GeV 


132 


2 500 


5.3 % 




2.6 


30 GeV 


75 


1 500 


5.0 % 




1.9 


40 GeV 


30 


400 


7.5 % 




1.5 



23 



450 




ID 


206 


400 








Mean 


109,0 


350 








RMS 


34,66 








ALLCHAN 


1263. 


300 








x'/ndf 13,28 


/ 2 


250 








Constant 


420,0 


200 








Mean 


101,1 








Sigma 


12,67 


150 










100 










50 


1 _LJ \ ^ 
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Figure 8: Distribution of the reconstructed m].jf, m,ass for the two best combinations 
in events containing four b-jets selected by algorithm B and for an integrated lumi- 
nosity of 3 • 10^ pb^^ (ATLAS). The top distribution is for H — > hh ^ bbbb signal 
events with hih = 300 GeV and tan/3 = 1.5, corresponding to iiih = 98 GeV. The 
bottom distribution is for the background from multi-jet production (solid), where 
the contributions from direct bb (cross-hatched) and bb + gb production (dashed) 
are also shown. 
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Figure 9: Same as Fig. for the m^bbb spectrum after applying a constraint, 
iribb = nih, on the two best bb combinations. 
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Figure 10: Distributions ofp!^, p^^^^, prp^ and prp* for signal (left) and background 
(right) after applying the selection B with a jet Yjt^ -threshold of 4.0 GeV, for an 
integrated luminosity of 3 ■ 10^ pb^^ (ATLAS). 
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Figure 11: Distributions ofp!^, p^^^^, prp^ and prp* for signal (left) and background 
(right) after applying the selection B with a jet Yjt^ -threshold of 20 GeV, for an 
integrated luminosity of 3 ■ 10^ pb^^ (ATLAS). 
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4 Conclusions 



This first study of event topologies containing four b-jets in the search for MSSM 
Higgs-boson decays has shown that it will be very difficult to cleanly extract the 
signal above the background over a significant region of the MSSM parameter space. 
The results of this note are much less optimistic than those from P], because the 
study reported in [§] did not consider all processes leading to multiple b-jets in the 
final state, nor did it consider the significant fraction of fake b-jets to be expected 
in such complex topologies. The study reported here was particularly difficult to 
bring to a conclusion since more than 2 ■ 10'' background events had to be generated 
and processed through ATLFAST. 

The bbH, bbA with H, A — > bb channel was considered for large values of mn, mA 
and tan (3. The strongly enhanced production cross-section of 3.5 pb for mH=500 GeV 
and tan/3 = 30 yields 10^ signal events for an integrated luminosity of 3 ■ 10^ pb~^. 
This large rate is however strongly reduced by the selection cuts and the b-tagging 
efficiencies, leading to roughly 5000 observable events in a mass bin of ±80 GeV 
around mn, mA. The QCD multi-jet background is initially several orders of mag- 
nitude higher and after the selection cuts and the b-tagging procedure a signal-to- 
background ratio of only a few percent can be achieved. The sensitivity to this 
channel is therefore rather weak. 

The H — > hh — > bbbb channel was studied for low tan (3 and mn = 300 GeV. 
Fully reconstructed h — >■ bb and H ^ hh decays in relatively narrow mass windows 
would give firm evidence for both h and H Higgs boson production. But in the 
presence of the huge QCD background, the sensitivity to this channel is also weak, 
weaker than that of the H ^ hh ^ bb77 channel. 

In both cases, at most only a fraction of the MSSM parameter space, which is al- 
ready covered by other decay modes (bbH, bbA, with H, A — >■ rr and H — hh — > bb77) 
is accessible to these channels containing four b-jets in the final state. Nevertheless, 
their observation would help in constraining the couplings and branching ratios of 
the MSSM Higgs bosons. 
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